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We study the thermodynamics of the Nambu-Jona-Lasinio model with Polyakov loops, 
where spontaneous chiral symmetry breaking and confinement are taken into account. We 
focus on the phase structure of the model and explore the susceptibilities associated with 
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corresponding order parameters under the mean-field approximation. 

§1. Introduction 



| Dynamical chiral symmetry breaking and confinement are related to global sym- 

metries of the QCD Lagrangian. While the chiral symmetry is exact in the limit of 
massless quarks the Z(3) center symmetry, which governs confinement, is exact in 
the opposite limit, i.e., for infinitely heavy quarks. The calculation in lattice gauge 
theory (LGT) with dynamical quarks may indicate that the global symmetries of the 
Lagrangian are still relevant even in presence of their explicit breaking. x ) Thus, to 
account for deconfinement in the presence of dynamical quarks, one can construct an 
effective Lagrangian in terms of the order parameters, the Ginzburg-Landau effective 
theory. 

Recently, gluon degrees of freedom were introduced in the Nambu-Jona-Lasinio 
■ (NJL) Lagrangian 2 ) through an effective gluon potential expressed in terms of Polyakov 

loops (PNJL model). 3 )' 4 ) The model has a non-vanishing coupling of the constituent 
quarks to the Polyakov loop and mimics confinement in the sense that only three- 
s' ' quark states contribute to the thermodynamics in the low-temperature phase. Hence, 
the PNJL model yields a better description of thermodynamic quantities near the 
phase transition than the NJL model. The PNJL model has been applied to explore 
the thermodynamical observables. 4 ) -10 ) Furthermore, due to the symmetries of the 
Lagrangian, the model belongs to the same universality class as that expected for 
QCD. Thus, the PNJL model can be considered as a testing ground for the critical 
phenomena related to the breaking of the global Z(3) and chiral symmetries. 

An extension of the NJL Lagrangian by coupling the quarks to a uniform tem- 
poral background gauge field, which manifests itself entirely in the Polyakov loop, 
has been proposed to account for interactions with the color gauge field in effective 
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chiral models , 3 )' 4 ) The PNJL Lagrangian is given by 

C = L N j L ^,$[A},§[A] ] T,^-U^[A],3[A]-,T). (1-1) 

The interaction between the effective gluon field and the quarks in the PNJL La- 
grangian is implemented by means of a covariant derivative in £njl- The effective 
potential IA of the gluon field in (|1T|) is expressed in terms of the traced Polyakov 
loop ^ and its conjugate ^ 

* = -TT^cL , $ = ^Tr c L f , (1-2) 

where L is a matrix in color space related to the gauge field by 

r/3 



L{x) =V exp 



/ drA^x, t) 
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(1-3) 



with V being the path (Euclidean time) ordering, and (3 = 1/T with A4 = iAo. In 
the heavy quark mass limit, the thermal expectation value of the Polyakov loop (#) 
acts as an order parameter of the Z(3) symmetry. 

The effective potential U(<P, 3) of the gluon field is expressed in terms of the 
Polyakov loops so as to preserve the Z(3) symmetry of the pure gauge theory. n ) We 
adopt an effective potential of the following form: 4 ) 

m^n = jm^.h^ +¥)+ h m ^ ( , 4) 

where the coefficients 6j are fixed by requiring that the equation of state obtained in 
pure gauge theory on the lattice is reproduced. 



§2. Susceptibilities and the phase structure 



In Fig. [U (left) W6 show th.6 ch-iral Xmm 

and Polyakov loop Xu = i®®} ~ (*) (#) 
susceptibilities computed at /U = in the PNJL model in the chiral limit within the 
mean field approximation. The chiral susceptibility exhibits a very narrow divergent 
peak at the chiral critical temperature T c h, while the peak of Xu is much broader 
and the susceptibility remains finite at all temperatures. This is due to the explicit 
breaking of the Z(3) symmetry by the presence of quark fields in the PNJL La- 
grangian. Nevertheless, Xu still exhibits a peak structure that can be considered as 
a remnant for the deconfinement transition in this model. Thus, for the parameters 
used in the model, the deconfinement transition, signaled by the broad peak in Xu-i 
sets in earlier than the chiral transition at vanishing net quark density. 

The peak positions of the Xmm and Xu susceptibilities can be used to determine 
the phase boundaries in the (T, //)-plane. In Fig. [T] (right) we show the resulting 
phase diagram for the PNJL model. The boundary lines of deconfinement and 
chiral symmetry restoration do not coincide. There is only one common point in 
the phase diagram where the two transitions appear simultaneously. Recent LGT 
results for two flavor both at vanishing and at finite quark chemical potential show 
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Fig. 1. Left: the chiral \ m m (dashed-line) and the Polyakov loop \u (solid- line) susceptibilities 
in the chiral limit as functions of temperature T for fi — 0. Right: the phase diagram of the 
PNJL model in the chiral limit. The solid (dashed) line denotes the chiral (deconfinement) 
phase transition respectively. The TCP (bold-point) is located at (T c = 157, (M c = 266) MeV. 



that deconfinement and chiral symmetry restoration appear along the same critical 
line. 1 ) In general it is possible to choose the PNJL model parameters such that the 
both critical temperatures coincide at fx = 0. 

From the phase diagram shown in Fig. [I] one finds that the slope of T^ ec as a 
function of n is almost flat, indicating that at low temperature the chiral phase tran- 
sition should appear much earlier than deconfinement. However, there are general 
arguments, that the deconfinement transition should precede restoration of chiral 
symmetry. 12 ^ In view of this, it seems unlikely that at T ~ the chiral symmetry 
sets in at the lower baryon density than deconfinement. In the PNJL model, the 
parameters of the effective gluon potential were fixed by fitting quenched LGT cal- 
culations. Consequently, the parameters are taken as independent on \i. However, it 
is conceivable that the effect of dynamical quarks can modify the coefficients of this 
potential, thus resulting in //-dependence of the parameters. The slope of Td ec as a 
function of \i could be steeper 0- Consequently, the effective Polyakov loop potential 
(|l-4j) with ^-independent coefficients should be considered as a good approximation 
only for fi/T < 1. 

While the susceptibilities Xmm and Xu exhibit expected behaviors associated 
with the phase transitions, the diagonal Polyakov loop susceptibilities, 

XB = <*">- <*> 2 , Xu = ($ 2 ) ~ ($) 2 , (2d) 

are negative in a broad temperature range above T c h. 14 ) This is in disagreement with 
recent lattice results, where in the presence of dynamical quarks xu 1S found to be 
always positive. ^ A possible reason for this behavior could be traced to the form of 
the effective Polyakov loop potential used in the Eq. (|1-4|) . 

Recently an improved effective potential with temperature-dependent coeffi- 
cients has been suggested 8 ^ in which the Haar measure in the group integral is ac- 
counted for. The improved potential indeed yields positive values for all the Polyakov 
loop susceptibilities. We note that the phase diagram calculated with the improved 
potential is similar to that obtained with the previous choice of the Polyakov loop 

*' Such a modification was explored 13 ' where explicit fx- and Nf - dependence in 6; is extracted 
from the running coupling constant a s , using arguments based on the renormalization group. 
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interactions, shown in Fig. [TJ 



§3. Conclusions 



We introduced susceptibilities related to the three different order parameters in 
the PNJL model, and analyzed their properties and their behavior near the phase 
transitions. In particular, for the quark-antiquark and chiral density-density corre- 
lations we have discussed the interplay between the restoration of chiral symmetry 
and deconfinement. We observed that a coincidence of the deconfinement and chiral 
symmetry restoration is accidental in this model. 

We found that, within the mean field approximation and with the polynomial 
form of an effective gluon potential the correlations of the Polyakov loops in the 
quark-quark channel exhibit an unphysical behavior, being negative in a broad 
parameter range. This behavior was traced back to the parameterization of the 
Polyakov loop potential. We argued that the Z(iV)-invariance of this potential and 
the fit to lattice thermodynamics in the pure gluon sector is not sufficient to pro- 
vide a correct description of the Polyakov loop fluctuations. Actually it was pointed 
out 8 ) that the polynomial form used in this work does not possess the complete 
group structure of color SU(3) symmetry. The improved potential yields a positive, 
i.e. physical, xu susceptibility, in qualitative agreement with LGT results. 
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